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Abstract. Based on a Riemann theta function and Hirota's bilinear form, a lucid 
' and straightforward way is presented to explicitly construct double periodic wave 

solutions for both nonlinear differential and difference equations. Once such a 
equation is written in a bilinear form, its periodic wave solutions can be directly 
obtained by using an unified theta function formula. The relations between the 
periodic wave solutions and soliton solutions are rigorously established. The 
efficiency of our proposed method can be demonstrated on a class variety of non- 
linear equations such as those considered in this paper, shall water wave equation, 
■ (2+l)-dimensional Bogoyavlenskii-Schiff equation and differential-difference KdV 
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1. Introduction 

The bilinear derivative method developed by Hirota is a powerful and direct approach 
to construct exact solution of nonlinear equations. Once a nonlinear equation is written 
in bilinear forms by a dependent variable transformation, then multi-soliton solutions are 
usually obtained [TJ-[6]. It was based on Hirota forms that Nakamura proposed a convenient 
way to construct a kind of quasi-periodic solutions of nonlinear equations [TJ |8], where the 
periodic wave solutions of the KdV equation and the Boussinesq equation were obtained. 
Such a method indeed exhibits some advantages. For example, it does not need any Lax 
pairs and Riemann surface for the considered equation, allows the explicit construction of 
multi-periodic wave solutions, only relies on the existence of the Hirota's bilinear form, as well 
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as all parameters appearing in Riemann matrix are arbitrary. Recently, further development 
was made to investigate the discrete Toda lattice, (2+l)-dimensional Kadomtsev-Petviashvili 
equation and Bogoyavlenskii's breaking soliton equation [9j-[14j. However, where repetitive 
recursion and computation must be preformed for each equation [7]-|14j. 

The motivation of this paper is to considerably improves the key steps of the above 
existing methods. To achieve this aim, we devise a theta function bilinear formula, which 
actually provides us a lucid and straightforward way for applying in a class of nonlinear 
equations. Once a nonlinear equation is written in bilinear forms, then the double periodic 
wave solutions of the nonlinear equation can be obtained directly by using the formula. 
Moreover, we propose a simple and effective method to analyze asymptotic properties of 
the periodic solutions. As illustrative examples, we shall construct double periodic wave 
solutions to the shall water wave equation, (2+l)-dimensional Bogoyavlenskii-Schiff equation 
and differential-difference KdV equation. 

The organization of this paper is as follows. In section 2, we briefly introduce a Hirota 
bilinear operator and a Riemann theta function. In particular, we provide a key formula 
for constructing double periodic wave solutions for both differential and difference equa- 
tions. As applications of our method, in sections 3-5, we construct double periodic wave 
solutions to the shall water wave equation, (2+l)-dimensional Bogoyavlenskii-Schiff equation 
and differential-difference KdV equation, respectively. In addition, it is rigorously shown that 
the double periodic wave solutions tend to the soliton solutions under small amplitude limits. 

2. Hirota bilinear operator and Riemann theta function 

To fix the notations we recall briefly some notions that will be used in this paper. The 
Hirota bilinear operators D x ,Dt and D n are defined as follows: 



D™D*f(x,t)-g(x,t) = (d x -d x ,) m {d t -d t ,) k f{x,t)g(x\t')\ xl -_ 
e 5D "f(n) ■ g(n) = e^ 9 ^ f(n)g(n')\ n , =n = f(n + 8)g(n - 6), 




Proposition 1. The Hirota bilinear operators D x ,Dt and D n have properties pQ-[6] 



cosh(5D n )e €l • e ?2 = cosh[5(z/i - i/ 2 )]e 6+6 , 
sinh^Aje* 1 • e 6 = sinh[5(i/i - ^ 2 )]e ?1+6 , 

where £j = ctyx + Wji + i/j-n + dj, and ay, Wj, f,-, cr,-, j = 1,2 are parameters and n € Z is a 
discrete variable. More generally, we have 



F(D x ,D t ,D n )e^ -e^ =F(a 1 -a 2 ,0J 1 -0J2,e^[5(u l -u 2 )})e^, 



(2.1) 



where F(D x ,D t ,D n ) is a polynomial about operators D x ,Dt and -D n . This properties are 
useful in deriving Hirota's bilinear form and constructing periodic wave solutions of nonlinear 
equations. 

In the following, we introduce a general Riemann theta function and discuss its periodic- 
ity, which plays a central role in the construction of periodic solutions of nonlinear equations. 
The Riemann theta function reads 



(£, r) = exp{27ri(£ + e)(m + s) — 7rr(m + s) 2 }. 



(2.2) 



Here the integer value m € Z, complex parameter s, e € C, and complex phase variables 
£ £ C; The r > which is called the period matrix of the Riemann theta function. 

In the definition of the theta function (2.2), for the case s = e = 0, hereafter we use 



0(£,r)=0 



(M=0(£ + e,r) 



(£, t) for simplicity. Moreover, we have 
Definition 1. A function g(t) on C is said to be quasi-periodic in t with fundamental 
periods Xi, • • ■ , G C, if Xi, • • • , Tfe are linearly dependent over Z and there exists a function 
G(yi, ■■■ ,yk), such that 

G(yi, ■■■ ,yj + Tj,--- ,y k ) = G(y x ,-- ■ ,y fc ), for all % € C, j = 1,- • • ,k. 

G{t,--- ,V ,t)=g{t). 

In particular, is called double periodic as k = 2, and it becomes periodic with T if and 
only if Tj = rrijT, j = 1, • • • , A;. □ 

Let's first see the periodicity of the theta function r). 

Proposition 2. [15J The theta function i?(£,t) has the periodic properties 

+ 1 + ir, t) = exp(-2vr^ + 7rr)tf (£, r). (2.3) 



We regard the vectors 1 and iTj as periods of the theta function r) with multipliers 1 
and exp(— 27ri£ + ttt), respectively. Here, it is not a period of theta function r), but it 
is the period of the functions <9| ln0(£, r), ln[#(£ + e, t)/0(£ + h, r)] and 0(£ + e, r)0(£ - 
e,r)/^ + ^r) 2 . 

Proposition 3. The meromorphic functions /(£) on C are as follow 
(«) /«) = 8flni»«,r), {€C, 

then in all three cases (i)-(iii), it holds that 

/(^+i + ir) = /(o, eec, 

that is, /(£) is a double periodic function with 1 and it. 
Proof. By using (2.3), it is easy to see that 



(2.4) 



tf(£ + l + iT,r) — ' 0(£, r ) ' 
or equivalently 

d 5 ln??(£ + l + iT,T) = -2m + d(]n0(Z,T). (2.5) 

Differentiating (2.5) with respective to £ again immediately proves the formula (2.4) for the 
case (i). The formula (2.4) can be proved for the cases (ii) and (iii) in a similar manner. □ 



Theorem 1. Suppose that ■& 
in which £ = ax + cut + vn + a. 



(£,r)andtf 



(£, r) are two Riemann theta functions, 



'hen Hirota bilinear operators D x ,Dt and D n exhibit the 



following perfect properties when they act on a pair of theta functions 



(&r)-0 



2 ^ 

M=0,1 



-H/2 



exp(5D n )$ 



(£,r)-0 



£ — £ 

-Ai/2 



(C,r) 

(2e,2r)| 5=0 
(2£,2r)| f=0 



e' + £ 
M/2 



(2.6) 



(2£,2r), 



£' + £ 

M/2 



(2.7) 



(2£,2r), 



exp(<5D n )tf 

ju=0,l 

where the notation X^=o l represents two different transformations corresponding to fx = 0, 1. 
The bilinear formula for t is the same as (2.6) by replacing d x with <9 t . 



In general, for a polynomial operator F(D x ,D t ,D n ) with respect to D x ,D t and D n , we 
have the following useful formula 



F(D x ,D tl D n )H 



(&r)-0 



1? 



V + £ 

M/2 



(2£,2t), (2.8) 



in which, explicitly 

C(e,e f ,n)= ^ F(7W)exp [-2vrr(m-^/2) 2 -27ri(m-^/2)(e' -e)] . ^ 9) 



mat 



where we denote vector M. = (47ri(m — [i/2)a, 4iri(m — fj,/2)u, exp[4-7ri(m — /u/2)<5f] ) . 
Proof. Making use of Proposition 1, we obtain the relation 



D x d 



(£,r)-0 



for) 



= D x exp{27rim'(^ + e') — 7rm /2 r} • exp{27rim(£ + e) — 7rm 2 r}, 

= 2-7ria(m' — m) exp \2m\m! + m)£ — 2m(m'e' + me) — 7rr[m /2 + m 2 ]} 

By shifting sum index as m = V — ml ', then 

A = 1ma(2m' - I') exp {2iril'£ - 2m[m'e' + (/' - m')e] - iTT[m' 2 + (/' - m') 2 ]} 

i'=a +A1 ^ 2 7ria(2m' - 2/ - n) exp{47ri£(Z + n/2) 

M=0,1 i,m'eZ 

-2vri[m / e / - (m - 21 - fi)e] - ir{m' 2 + (m 1 -21- /u) 2 ]r} 
Finally letting m! = k + /, we conclude that 



47ria[fc - /ii/2] exp{-27ri(fc - /i/2)(e' - e) - 2nT{k - n/2) 2 } 



x ^ exp{27ri(Z + /Li/2)(2£ + e' + e) - 27rr(Z + /i/2) 2 



U*=o,i 



-/i/2 



d 



e' + e 
/i/2 



(2£,2t), 



(2^,2r)| £=0 

by using the following relations 

k + l = (k - ft/2) + + /x/2), k-l-n=(k- /i/2) - (I + /i/2). 

In a similar way, we can prove the formula (2.7). The formula (2.8) follows from (2.6) 
and (2.7). □ 

Remark 1. The formulae (2.8) and (2.9) show that if the following equations are satisfied 

C(e,e', f i) = 0, (2.10) 
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for jjL = 0, 1, then $ 
equation 



(f,r) and ■d 



(£,t) are periodic wave solutions of the bilinear 



F(D x ,D t ,D n )$ 



The formula (2.10) contains two equations which are called constraint equations. This formula 
actually provides us an unified approach to construct double periodic wave solutions for both 
differential and difference equations. Once a equation is written bilinear forms, then its 
periodic wave solutions can be directly obtained by solving system (2.10). 

Theorem 2. Let C(e, s',fJ>) and F(D X , Dt, D n ) be given in Theorem 1, and make a 
choice such that e' — e = ±1/2. Then 

(i) If F(D X , D t , D n ) is an even function in the form 

F(-D x , -D t , -D n ) = F(D X , D t , D n ), 
then C(e,s',fj,) vanishes automatically for the case \i = 1, namely 

C(£,e',n) = 0, for fj, = 1. (2.11) 

(ii) If F(D X , Dt, D n ) is an odd function in the form 

F(-D x , -D t , -D n ) = -F(D X , D t , D n ), 

then C(e,£',/x) vanishes automatically for the case \i = 0, namely 

C{e,e',v) = 0, for^ = 0. (2.12) 

Proof. We are going to consider the case where F(D X , Dt, D n ) is an even function and 
prove the formula (2.11). The formula (2.12) is analogous. Making transformation m = 
—fh + fx, and noting F(D X , D t , D n ) is even, we then deduce that 

C(e, e', n) = F(-M) exp [-27rr(m - fi/2) 2 + 2m{fh - n/2){e' - e)] 

fh&L 

= C(e, e' , n) exp [4iri(fh — n/2)(e' — e)] 

= C(e, e' , fx) exp (±2nifh) exp (±7ri/i) = —C(e, e' , fx), 

which proves the formula (2.11). □ 

Corollary 1. Let e'j — Ej = ±1/2, j = l, -- ,N. Assume F{D x ,D t ,D n ) is a linear 

combination of even and odd functions 



F(D X , D tl D n ) = Fi(D x , D t , D n ) + F 2 (D X , D t , D n ), 
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where Fi(D x ,D t ,D n) is even and F 2 (D X , Dt, D n ) is odd. In addition, C(£,£',/i) correspond- 
ing (2.9) is given by 

C(s,e',n) = Ci(e,e',/i) + C 2 (e,e',)u), 

where 

Ci(e,e', / u)= ^ Fi(7W)exp[-27rT(m-/i/2) 2 -27ri(m-^/2)(e'-e)] , 
C 2 (e, £',//) = F 2 {M)exp[-2TTT(m- fi/2) 2 -2TTi(m- fj,/2)(e' - e)] . 

Then 

C(e,e',/i) = C7 2 (e,e', m) for ^ = 1, (2.13) 

C(e,6/,n) = Ci(e,e?,n), for /x = 0. (2.14) 

Proof. In a similar to the proof of Theorem 2, shifting sum index as m = —in + and 
using Fi(D x , D t , D n ) even and F 2 (D X , D t , D n ) odd, we have 
C{e,e',n) = Ci(£,e',n) + C 2 (e,e',/i) 

(2.15) 

= [Ci(E,e',n) - C 2 (e,£ , ,n)] exp(±mpi) . 
Then for \x = 1, the equation (2.15) gives 

Ci(e,e',At) = 0, 

which implies the formula (2.13). The formula (2.14) is is analogous. □ 

The theorem 2 and corollary 1 are very useful to deal with coupled Hirota's bilinear 
equations, which will be seen in the following section 4. 

3. The shall water wave equation 

The shall water wave equation takes the form |16j 

/•oo 

u t - u xxt - 3uu t + 3u x u t dx + u x = 0, (3.1) 

J x 

which is like to the KdV equation in the family of shall water wave equations. Hirota and 
Satsuma obtained soliton solutions of the equation by means of bilinear method [T7j. Here 
we construct its a double periodic wave solution and show that the one-soliton solution can 
be obtained as limiting case of the double periodic solution. 

To apply the Hirota bilinear method for constructing double periodic wave solutions of 
the equation (3.1), we consider a variable transformation 

u = 2ti%hif(x,t). (3.2) 



Substituting (3.2) into (3.1) and integrating with respect to x, we then get the following 
Hirota's bilinear form 

F(D x ,D t )f ■ f = (D x D t + D 2 X - D t Dl + c)f ■ f = 0, (3.3) 

where c is an integration constant. In the special case of c = 0, starting from the bilinear 
equation (3.3), it is easy to find its one-soliton solution 

ui = 2^1n(l + e"), (3.4) 

with phase variable r\ = px + ^rzft + 7 for every p and 7. 

Next, we turn to see the periodicity of the solution (3.2), the function / is chosen to be 
a Riemann theta function, namely, 

f(x,t) = 0(£,T), (3.5) 

where phase variable £ = ax + ojt + a. With Proposition 3, we refer to 

u = 2d 2 x \nd{^,T) = 2a 2 df lntf(£,r), (3.6) 

which shows that the solution u is a double periodic function with two fundamental periods 
1 and it. 

We introduce the notations by 

A = e -W2 5 j A) = 0(2£ 2r ) = £ A 4m ' exp(4™0, 

m€Z 

(3.7) 

(2£, 2r) = ^ A^" 1 - 1 ) 2 exp[2m(2m - 

m£Z 



^ 2 (e,A) = T? 





-1/2 



where the phase variable £ = ax + + a. 

Substituting (3.5) into (3.3), using formula (2.10) and (3.7) leads to a linear system ( 
corresponding to \i = and fj, = 1, respectively) 

[<(0, A)a + ^ 4) (0, A)a 4 ]w + #i(0, A)c + <(0, A)a 2 = 0, 

(3.8) 

k? 2 '(0, A)a + 4 4) (°> A)a 4 ]w + 2 (O, A)c + tf 2 '(0, A)a 2 = 0, 
where we have denoted the derivative of A) at £ = by notations 

(Qj A) = ^M| H) i = 1, 2; /c = 1, 2, 3, 4. 
This system admits an explicit solution (oj,c). In this way, we obtain an explicit periodic 



wave solution (3.6) with parameters uj, c by (3.8), while other parameters a, a, r, a are free. 

In summary, double periodic wave (3.6) possesses the following features: (i) It is is one- 
dimensional, i.e. there is a single phase variable £. Moreover, it has two fundamental periods 



1 and it in phase variable £, but it need not to be periodic in x and t . (ii) It can be viewed 
as a parallel superposition of overlapping one-soliton waves, placed one period apart. 

In the following, we further consider asymptotic properties of the periodic wave solu- 
tion. Interestingly, the relation between the periodic wave solution (3.6) and the one-soliton 
solution (3.4) can be established as follows. 

Theorem 3. Suppose that the vector (oj, c) is a solution of the system (3.8), and for the 
periodic wave solution (3.6), we let 

p J + ttt 

a = ^V a = ^T' (3 ' 9) 

where the p and 7 are given in (3.4). Then we have the following asymptotic properties 

P 

c — > 0, 2irit; — txt — > ri = px H — » 1 + 7, 

p l — 1 

#(Z,t) — ►l + e'', as A^O. 
In other words, the double periodic solution (3.10) tends to the soliton solution (3.4) under 
a small amplitude limit, that is, 

u — >ui, as A -> 0. (3.10) 

Proof. Here we will directly use the system (3.8) to analyze asymptotic properties of the 
periodic solution (3.6). Since the coefficients of system (3.8) are power series about A, its 
solution (u, c) also should be a series about A. We explicitly expand the coefficients of system 
(3.8) as follows 

0i (0, A) = 1 + 2A 4 + • • • , <(0, A) = -32tt 2 A 4 + • • • , 

tfJ 4) (0,A) = 512tt 4 A 4 + • • • , #2(0, A) = 2A + 2A 9 + ••• (3.11) 

0£(O, A) = -8vr 2 A + • • • , 4 4) (0, A) = 32tt 4 A + • • • . 

Let the solution of the system (3.8) be in the form 

uj = u) + u>\\ + u> 2 X 2 -\ = ojo + o(A), 

2 (3-12) 
c = c + ciA + c 2 A H = c + o(A). 

Substituting the expansions (3.11) and (3.12) into the system (3.8) (the second equation 

is divided by A ) and letting A — > 0, we immediately obtain the following relations 

Co = 0, (-87r 2 a + 32vr 4 a 3 )w - 8vr 2 a 2 = 0, 

which implies 
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Combining (3.12) and (3.13) then yields 



c — > 0, 2iriuj — > ln _ m x9 = — -, as A — > 0. 



2-nia. p 
(2Tria) 2 — 1 p 
Hence we conclude 

£ = 27TZ£ — TTT = px + 2Tliujt + 7 



V 

px + -5 -t + 7 = 77, as A ->■ 0. 

p z — 1 



(3.14) 



It remains to consider asymptotic properties of the periodic wave solution (3.6) under the 
limit A — > 0. By expanding the Riemann theta function z?(£,t) and using (3.14), it follows 
that 

0(£ r) = 1 + A 2 (e 2 ^ + e~ 2 ^) + A 8 (e 4 ^ + e" 4 ^) + • • • 

= i + J + A 4 (e"« + e 2 «) + A 12 (e" 2 « + e 3 «) + • • • 

— ^1 + e 5 — >l + e v , as A ->■ 0, 
which together with (3.6) leads to (3.10). Therefore we conclude that the double periodic 

solution (3.6) just goes to the one-soliton solution (3.4) as the amplitude A — > 0. □ 
4. The modified Bogoyavlenskii-Schiff equation 

We consider (2+l)-dimensional modified Bogoyavlenskii-Schiff equation [18] 

u t - iu 2 u z - 2u x d~ l (u 2 ) z + u xxz = 0, (4.1) 

which was deduced from the Miura transformation }19| . Equation (4.1) is reduced to the 
modified KdV equation in the case of x = z. 

We shall construct a double periodic wave solution to the equation (4.1) by using Theorem 

1 and 2. The equation (4.1) can be described by a coupled system 

u = ip x , 

Pxx + ^l + c = 0, (4.2) 

ipt + 2ip x p xz + ip z {p xx + ip 2 x + c) + tp xxz = 0. 
We perform the dependent variable transformations 

u = ip x = d x In ( - J , p = ln(fg), (4.3) 



(4.4) 



JJ, 

then equation (4.2) is reduced to the following bilinear form 
F(D x )f ■ g = (D 2 X + c)f ■ g = 0, 

G(D t , D x ,D z )f ■g = (D t + D 2 X D Z + cD z )f ■ g = 0, 
where c is a constant. The equation (4.4) is a type of coupled bilinear equations, which is 
more difficult to be dealt with than the single bilinear equation (3.3) due to appearance of 
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two functions and two equations. We will take full advantages of Theorem 2 to reduce the 
number of constraint equations. 

Now we take into account the periodicity of the solution (4.3), in which we take / and g 

as 

/ = 0(£ + e,r), g = #(Z + h,T), e,heC, (4.5) 

where phase variable £ = ax + (3z + cot + a. By means of Proposition 3, we find that the 
solution 

#(g + e,T) 
u = adt In — — — - 

6 #(Z + h,T) 

is a double periodic function with two fundamental periods 1 and it. 

In the special case of c = 0, the equation (4.2) admits one-soliton solution 

1 + e v 

i*i = Or In— — , (4.6) 

where rj = px + qy — p 2 qt + 7 for every p, q and 7. 
We take e = 0, h = 1/2 in (4.5), and therefore 

/ = t) = ^2 exp(27rin£ - TrmV), 



1/2] 




(£> r ) = S ex P( 2 ™™(£ + I/ 2 ) - vrmV) (4.7) 
= X] ( _1 ) m ex P( 2 ™K - 7rm 2 r). 

mgZ 

Due to the fact that F{D X ) is an even function, its constraint equations in the formula 
(2.10) vanish automatically for fi = 1. Similarly the constraint equations associated with 
G(D t , D x , D z ) also vanish automatically for ji = 0. Therefore, the Riemann theta function 
(4.6) is a solution of the bilinear equation (4.4), provided the following equations 

<(0,A)a 2 +t?i(0,A)c = 0, 

(4.8) 

tf' 2 (0, A)w + 0^(0, A)/3c + <(0, A)a 2 /3 = 0, 
where we introduce the notations by 

A = e- 7 ^ 2 , A) = 0(2£, 2r) = £ A 4 ™ 2 exp(4mm0, 



^(e,A) = 1 ? 



1/2" 



-1/2 



(2£,2r) = ^(-lf^ 2 ™- 1 ) exp[2OT(2m - 

It is obvious that equation (4.8) admits an explicit solution uj and c. In this way, a periodic 
wave solution reads 
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where parameters to and c are given by (4.11), while other parameters a,/3,T, a are free. In 
summary, double periodic wave (4.9) has the following features: (i) It is one-dimensional and 
has two fundamental periods 1 and it in phase variable £. (ii) It can be viewed as a parallel 
superposition of overlapping one-soliton waves, placed one period apart. 

In the following, we further consider asymptotic properties of the double periodic wave 
solution. The relation between the periodic wave solution (4.9) and the one-soliton solution 
(4.6) can be established as follows. 

Theorem 4. Suppose that the vector (w,c) T is a solution of the system (4.8). In the 
periodic wave solution (4.9), we choose parameters as 

a - — f3 - — a - 7 + — (4 10) 

where the p, q and 7 are the same as those in (4.6). Then we have the following asymptotic 
properties 

c^O, ^1±^, /—n + e", g^\- e \ as A -> 0. 
2m 

In other words, the double periodic solution (4.9) tends to the one-soliton solution (4.6) under 
a small amplitude limit , that is, 

u — > ui, as A -> 0. (4-11) 

Proof. Here we will directly use the system (4.8) to analyze asymptotic properties of 
periodic solution (4.9). We explicitly expand the coefficients of system (4.8) as follows 
0i (0, A) = 1 + 2A 4 + • • • , <(0, A) = -32vr 2 A 4 + • • • , 

(4-12) 

tf' 2 (0, A) = —AiriX + UniX 9 + • • • , i?' 2 "(0, A) = 16yr 3 iA - 487r 3 iA 9 + • • • , 
Suppose that the solution of the system (4.8) is of the form 

oj = u) + ujiX + ijj 2 X 2 H =^o + o(A), 

(4-13) 

c = c + ciA + c 2 A H = c + o(A). 

Substituting the expansions (4.12) and (4.13) into the system (4.8) and letting A — > 0, 
we immediately obtain the following relations 

c = 0, -4ttiu; + 16ir 3 ia 2 (3 = 0, 

which has a solution 

c = 0, w = 47rV/3. (4.14) 
Combining (4.13) and (4.14) leads to 

c — > 0, 2ttiuj — > 8ir 3 ia 2 f3 = -p 2 q, as A -> 0, 
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or equivalently 

£ = 2-7ri£ — ivT = px + qy + 2iriujt + 7 

(4.15) 

— > px + qy — p qt + j = rj, as A — > 0. 
It remains to identify that the periodic wave (4.9) possesses the same form with the 

one-soliton solution (4.6) under the limit A — > 0. For this purpose, we start to expand the 

functions / and g in the form 

/ = 1 + A 2 (e 2 ^ + e~ 2 ^) + A 8 (e 47ri? + e" 4 ^) + • ■ ■ . 
g = l _ A 2 (e 2 ™ f + e~ 2 ^) + A 8 (e 4 ^ + e" 4 ^) + • • • . 
By using (4.13)-(4.15), it follows that 

/ = l + J + A 4 (e"« + e 2 «) + A 12 (e" 2 « + e 3 «) + • • • 

— ^l + e ? " — >l + e v , as A ->• 0; 

(4.16) 

g = 1 _ e i + A 4 (e 2 « - e"*) + A 12 (e~ 2 « - e 3 «) + • • • 

— ► 1 - — ► 1 - e v , as A -> 0. 
The expression (4.11) follows from (4.16), and thus we conclude that the double periodic 
solution (4.9) just goes to the one-soliton solution (4.6) as the amplitude A — > 0. □ 

5. The differential-difference KdV equation 

We consider differential-difference KdV equation 
d ( u(n) 



«i(") = - ( 5 - 2 ) 



u(„-l/2)-u(n+l/2). (5.1) 

at \ 1 + u(n) 

Hirota and Hu have found its soliton solutions and rational solutions [201 121| , among them 
one-soliton solution reads 

(l + e r?+p/2)( 1 + e r,-p/2) 

where rj = pn — smh(p/2)t + 7 for every p and 7. 

We shall construct a periodic wave solutions to the equation (5.1) by using Theorem 1. 
By means of a variable transformation 

f{n + l/2)/(n - 1/2) 
/(») 

the equation (5.1) is reduced to the bilinear equation 



u(n) = jT^Tg 1, (5.3) 



1 11 

sinh(- A) A + 2 sinh(-A) sinh(-A) + c 



/(n) ■ f[n) = 0, (5.4) 



where c is a constant. 
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Now we take into account the periodicity of the solution (5.3), in which we take f(n) = 
r), where phase variable £ = vn + ujt + a. Then solution (5.3) is written as 

"(0 = «(«) = ^-^1 ( 5 - 5 ) 

By means of Proposition 2, it is easy to deduce that u n is a double periodic function with 
two fundamental periods 1 and it. 

Substituting (5.5) into (5.4) and using formula (2.10) leads to a linear system 

sinh(- J D n )??' 1 (0, A)w + i?i(0, A)c + sinh( ^ D n ) sinh(-D„)i?i(0, A) = 0, 

1 1 1 ^ 

sinh(-D„)i9 2 (0, A)w + ?9 2 (0, A)c + sinh(- J D n ) sinh(-D n )i? 2 (0, A) = 0, 

where i?i(£, A) and $2(£, A) are the same as those in (3.7) with £ = z/n + ojt + <r. By using the 
solution u; and c of system (5.6), a periodic wave solution is obtained by (5.5). 

In the following, we further consider asymptotic properties of the double periodic wave 
solution. The relation between the periodic wave solution (5.5) and the one-soliton solution 
(5.2) can be established as follows. 

Theorem 5. Suppose that the vector (uj,c) t is a solution of the system (5.6). In the 
periodic wave solution (5.5), we choose parameters as 

v = ^v a = ^r> (5 - 7) 

where the p and 7 are the same as those in (5.2). Then we have the following asymptotic 
properties 

c^O, ^— >^±II, 0(£ jT )_n + e", as A -> 0. 

In other words, the periodic solution (5.5) tends to the one-soliton solution (5.2) under a 
small amplitude limit , that is, 

u(n) — >ui(n), as A — > 0. (5.8) 

Proof. Here we will directly use the system (5.6) to analyze asymptotic properties of 
periodic solution (5.5). We explicitly expand the coefficients of system (5.6) as follows 

i5i(0,A) = 1 + 2A 4 + --- , sinh(- J D n )i?' 1 (0,A) = 87risinh(ivri/)A 4 + • • • , 

sinh(-D n ) sinh(-D n )i?i(0, A) = 2 wah(iwi/) sinh(2i7rz/)A 4 + • • • , 

4 2 1 (5-9) 

d 2 (0, A) = 2A + 2A 9 + • • • , sinh(-D n ) 1 ? 2 (0, A) = Am sm\i(mv/2)\ + • • • , 

sinh(-D n ) sinh(-D n )i?2(0, A) = 2 sinh(i7r^) sinh(i7rz>/2)A + • • • . 
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Suppose that the solution of the system (5.6) is of the form 

oj = u>q+ ui\ + w 2 A 2 H = ojq + o(A), 

(5-10) 

c = co + ciA + C2A H = co + o(A). 

Substituting the expansions (5.9) and (5.10) into the system (5.6) and letting A — > 0, 
we immediately obtain the following relations 

Co = 0, 4-iri sinh(i7rz//2)cjo + 2 sinh(i7r^/2) sinh(i7rz/) = 0, 

which implies 

Co = 0, wq = — - — :sinh(i7r^). (5-H) 
Combining (5.9) and (5.10) leads to 

c — > 0, 2irioj — > — smh.(i-iru) = — sinh(p/2), as A — >■ 0, 

or equivalently 

£ = 2-irit; — ttt = pn + 2niujt + 7 

(5.12) 

— > pn — sinh(p/2)t + 7 = 77, as A — > 0. 
It remains to consider asymptotic properties of the periodic wave solution (5.5) under 

the limit A — >■ 0. By expanding the Riemann theta function i9(£,t), it follows that 

, r ) = 1 + J + A 4 (e-« + e 2 «) + A 12 (e" 2 « + e 3 «) + • • • 

— >• 1 + J- — >l + e v , as A ->• 0, 
which together with (5.5) lead to (5.8). Therefore we conclude that the periodic solution 
(5.5) just goes to the one-soliton solution (5.2) as the amplitude A — > 0. □ 
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